Cosmic Necklaces and Ultrahigh Energy Cosmic Rays 

Veniamin Berezinsky^, Alexander Vilenkin* 
^ INFN, Laboratori Nazionali del Gran Sasso, 67010 Assergi (AQ) Italy 
* Institute of Cosmology, Department of Physics and Astronomy , Tufts University, Medford, MA 

02155, USA 

' Cosmic necklaces are hybrid topological defects consisting of monopoles 



and strings, with two strings attached to each monopole. We argue that the 
cosmological evolution of necklaces may significantly differ from that of cosmic 



< 

• strings. The typical velocity of necklaces can be much smaller than the speed 

(N 

of light, and the characteristic scale of the network much smaller than the 

T 1 ■ 

> 

' horizon. We estimate the flux of high-energy protons produced by monopole 

CN ' annihilation in the decaying closed loops. For some reasonable values of the 

o 

0^ 



parameters it is comparable to the observed flux of ultrahigh-energy cosmic 
rays. 
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^ ' The observation of cosmic ray particles with energies higher than 10^^ GeV |l| gives 

^ i a serious challenge to the known mechanisms of acceleration. The shock acceleration in 
5t ! different astrophysical objects typically gives maximal energy of accelerated protons less than 
(1 — 3) ■ 10^° GeV [0]. The unipolar induction can provide the maximal energy 1 • 10^^ GeV 
only for the extreme values of the parameters Much attention has recently been given 
to acceleration by ultrarelativistic shocks @], [^. The particles here can gain a tremendous 
increase in energy, equal to F^, at a single reflection, where F is the Lorentz factor of the 
shock. However, it is known (see e.g. the simulation for pulsar relativistic wind in 0) that 
particles entering the shock region are captured there or at least have a small probability to 
escape. 

Topological defects (for a review see [0) can naturally produce particles of ultrahigh 
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energies (UHE) well in excess of those observed in cosmic rays (CR). In most cases the 
problem with topological defects is not the maximal energy, but the fluxes. 

Cosmic strings can produce particles when two segments of string come into close contact, 
as in cusp events . When the distance between two segments of the cusp becomes of the 
order of the string width, the cusp may "annihilate" turning into high energy particles. 
However, the resulting cosmic ray flux is far too small 

Superconducting strings [|3 appear to be much better suited for particle production. 



Moving through cosmic magnetic fields, such strings develop electric currents and copiously 
produce charged heavy particles when the current reaches certain critical value. The CR flux 
produced by superconducting strings is affected by some model-dependent string parameters 
and by the history and spatial distribution of cosmic magnetic fields. Models considered so 
far failed to account for the observed flux [|ll|,|l3 . 



Monopole-antimonopole pairs {MM) can form bound states and eventually annihilate 
into UHE particles ]13|], For an appropriate choice of the monopole density um-, this 
model is consistent with observations; however, the required (low) value of hm may be 
difficult to explain. 

We shall consider here another potential source of UHE CR, the topological defects 
which have not been much studied so far: cosmic necklaces. Such defects can be formed 
in a sequence of symmetry breaking phase transitions G ^ H x U{1) ^ H x Zi- If 
the group G is semisimple, then the first phase transition produces monopoles, and at the 
second phase transition each monopole gets attached to two strings, with its magnetic fiux 
channeled along the strings. The resulting necklaces resemble "ordinary" cosmic strings with 
monopoles playing the role of beads. "Realistic" particle physics models with necklaces can 



readily be constructed |15 



The evolution of necklaces is rather complicated, and its analysis would require high- 
resolution numerical simulations. Here we shall attempt only to indicate the relevant physical 
processes and to give very rough estimates for the efficiency of some of these processes. 

The monopole mass m and the string tension /i are determined by the corresponding 



symmetry breaking scales, rjs and r]m {Vm > Vs) , 



(1) 



Here, e is the gauge coupling. The mass per unit length of string is equal to its tension, 
Each string attached to a monopole pulls it with a force F — /lin the direction of the string. 
The monopole radius 5m and the string thickness Sg are typically of the order 5m ~ {^Vm)~^, 



Monopoles are formed at a temperature Tm ~ ?7m- Their initial average separation, d, 
can range from dm (for a second-order phase transition) to the horizon size (for a strongly 
first-order transition) . The monopoles are diluted by the expansion of the universe, so that d 
grows as d (X T~^. There is some additional decrease in the monopole density, and associated 
increase in d, due to MM annihilation. The latter process, however, is rather inefficient. 

At the second phase transition, each monopole gets attached to two strings, resulting 
in the formation of necklaces. There will be infinite necklaces having the shape of random 
walks and a distribution of closed loops. The two strings attached to a monopole are pulling 
it with an equal force; hence, there is no tendency for a monopole to be captured by the 
nearest antimonopole , unless their separation is comparable to the string thickness, Sg- 

An important quantity for the necklace evolution is the dimensionless ratio 



The average mass per unit length of necklaces is (r + l)/i. The initial value of r can be large 
(r S> 1) or small (r ^ 1), depending on the nature of the two phase transitions. 

We expect the necklaces to evolve in a scafing regime. If ^ is the characteristic length 
scale of the network, equal to the typical separation of long strings and to their characteristic 
curvature radius, then the force per unit length of string is / ~ /i/^, and the acceleration is 
a ~ (r -I- 1)~^^~^. We assume that ^ changes on a Hubble time scale ~ t. Then the typical 
distance travelled by long strings in time t should be ~ ^, so that the strings have enough 
time to intercommute in a Hubble time. This gives at^ ~ ^, or 



-1 



r = m/iid, 



(2) 
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e ~ (r + l)-'/H. (3) 

The typical string velocity is f ~ (r + 

For r <C 1 the monopoles are subdominant, and the string evolution is essentially the 
same as that of 'ordinary' strings without monopoles. The opposite case r ^ 1 is much 
different: the string motion is slow and their average separation is small. Like ordinary 
strings, cosmic necklaces can serve as seeds for structure formation. Significant quantitative 
changes in the corresponding scenario can be expected for r ^ 1. 

Disregarding MM annihilation, the evolution of r(t) can be analyzed using the energy 
balance equation 

E = ~PV - Eg. (4) 

Here, E is the energy of long necklaces in a co-moving volume V, P is the effective pressure, 
and Eg is the rate of energy loss by gravitational radiation from small-scale wiggles on 
long strings. If the scale of the wiggles is set by the gravitational back-reaction, then the 
strings radiate a substantial part of their energy in a Hubble time ll^Jl^] , and we can write 
Eg = KgNm/rt where N is the number of monopoles in volume V and ~ 1. The effect of 
loop formation is not relevant for the evolution of r{t) and has not been included in Eq.(^). 

For r <C 1, the effect of monopoles on the string dynamics is negligible, and we can 
write P = {N'm/3Vr){2v'^ — 1), where v is the rms string velocity. Then, with a power-law 
expansion a{t) cx t", we obtain the following equation for r(t): 

where Kg = z/(l — 2f^). The first term on the rhs of Eq.(|^) describes the string stretching 
due to expansion of the Universe while the second term describes the competing effect of 



string shrinking due to gravitational radiation |18]. In this regime, we can use the values 



of V from the string simulations [|1^]: v = 0.43 in the radiation era and v = 0.37 in the 
matter era. The corresponding values of Kg are, respectively, 0.07 and 0.14. Our estimate 
for Kg is Kg ~ 1, so it seems reasonable to assume that Kg > Ks. The solution of Eq.@ 



is r{t) oc i'^g-'^s^ suggesting that if r is initially small, it will grow at least until it reaches 
values r ~ 1. 

An equation similar to (pi) can also be written for r > 1, but in this case the results of 



numerical simulations |T9[ can no longer be used, and the relative magnitude of Ks and Kg 
cannot be assessed. Order-of- magnitude estimates suggest ~ k^, ~ 1, and in this paper 
we shall assume that Kg > Ks, so that r{t) is driven towards large values, r ^ 1. 

As r grows and monopoles get closer together, MM annihilation should become im- 
portant at some point. In any case, the growth of r should terminate at the value 
'^max ~ P-/rn6s ~ Vm/Vsj wheu the monopole separation is comparable to the string thickness 
6s. It is possible that annihilations will keep r at a much smaller value. For example, if 
monopoles develop appreciable relative velocities along the string, they may frequently run 
into one another and annihilate. The terminal value of r cannot be determined without 
numerical simulations of network evolution; here we shall assume that r ^ 1. 

Self-intersections of long necklaces result in copious production of closed loops. For r ^ 1 
the motion of loops is not periodic, so loop self-intersections should be frequent and their 
fragmentation into smaller loops very efficient. A loop of size i typically disintegrates on a 
timescale r ~ r~^/^£. All monopoles trapped in the loop must, of course, annihilate in the 
end. 

Annihilating MM pairs decay into Higgs and gauge bosons, which we shall refer to 
collectively as X-particles. The rate of X-particle production is easy to estimate if we note 
that infinite necklaces lose a substantial fraction of their length to closed loops in a Hubble 
time. The string length per unit volume is ~ and the monople rest energy released per 
unit volume per unit time is r/x/^^t. Hence, we can write 

dnx T'^/i 



(6) 



dt t^mx ' 

where mx ~ erjm is the X-particle mass and we have used Eq.(|^) . 

In the extreme case of r ~ Vmax ~ Vm/Vs, Eq-® gives the rate of X-particle production 
which does not depend on the string scale t]s. It is possible that the evolution of r{t) is 



actually saturated in this regime. 

X-particles emitted by annihilating monopoles decay into hadrons, photons and neutri- 
nos, which contribute to the spectrum of cosmic ultra-high energy radiations. In particular 
the diffuse flux of ultra-high energy protons can be evaluated as 

where dnx/dt is given by Eq.(||), Xp{E) is the attenuation lengh for ultra-high energy pro- 
tons due to their interaction with microwave photons and WN{mx,x) is the fragmentation 
function of X-particle into nucleons of energy E = xmx- 

The fragmentation function is calculated using the decay of X-paricle into QCD partons 
(quark, gluons and their supersymmetric partners) with the consequent development of the 
parton cascade. The cascade in this case is identical to one initiated by e^e^ -anihilation. 
We have used the fragmentation function in the gaussian form as obtained in MLLA approx- 
imation in ||20| and [^]. Additionally, we took into account the supersymmetric corrections 
to the coupling constant at large Q^. The details will be described elsewhere. Here we 
shall give only the explicit form of the fragmentation function we used: 

W^Ar(mx, x) = —exp ( 1 , (8) 



where 



6 V A 

X = E/mxi Xm = (A/mx)^''^, A = 0.234 GeV with the normalization constant to be 
found from energy conservation assuming that about 10% of initial energy {mx) is trans- 
ferred to nucleons. 

For attenuation length of UHE protons due to their interactions with microwave photons 
we used the calculations described in the book [||. 

Note that in our calculations the UHE proton flux is fully determined by only two 
parameters, r^/x and mx- The former is restricted by low energy diffuse gamma-radiation. 
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It results from e-m cascades initiated by high energy photons and electrons produced in the 
decays of X-particles. 

The cascade energy density predicted in our model is 



Ulr 



where to is the age of the Universe (here and below we use h = 0.75), z is the redshift and 
/,r ~ 1 is the fraction of energy transferred to pions. In Eq. (P) we took into account that half 
of the energy of pions is transferred to photons and electrons. The observational bound on the 
cascade density, for the kind of sources we are considering here, is ||2^ Ucas ~ 10~^ eV/cm^. 
This gives a bound on the parameter r^/i. 

In numerical calculations we used r^/x = 1 x 10^^ GeV^, which results in Ucas = 
5.6 • 10~^ eV/cm^, somewhat below the observational limit. Now we are left with one 
free parameter, mx, which we fix at 1 • 10^'' GeV. Note that with this value, the maximum 
energy of protons is not very high: E^ax ~ 10^'^ GeV. The calculated proton flux is pre- 
sented in Fig.l, together with a summary of observational data taken from ref. These 
data are usually interpreted as indicating the presence of a new component at energy higher 
than 1 • 10^° GeV. One cannot claim that our predicted flux gives a good fit to the data for 
this component, but the discrepancy does not exceed 2a. 

Let us now turn to the calculations of UHE gamma-ray flux from the decays of X- 
particles. The dominant channel is given by the decays of neutral pions. The flux can be 
readily calculated as 

A(^) = i^^^7(^)A^7(^), (10) 

where dnx/dt is given by Eq.(|^), X-yiE) is the absorption length of a photon with energy E 
due to e'^e' pair production on background radiation and Ny{E) is the number of photons 
with energy E produced per one decay of X-particle. The latter is given by 



The normalization constant K^^o is again found from the condition that neutral pions take 
away /7r/3 fraction of the total energy mx- 

An important point in our calculations was accounting for the absorption of UHE photons 
due to e+e~ production on background radiation. At energy E>1- 10^° GeV the dominant 
contribution to the absorption comes from the radio background. The significance of this 
process was first noticed in [|^(see also book |^). New calculations for this absorption were 



recently done . We have used the absorption lengths from this work. 



The calculated flux of gamma radiation is presented in Fig. 1 by the curve labelled 
7. One can see that at ~ 1 • 10^^ GeV the gamma ray flux is considerably lower than 
that of protons. This is mainly due to the difference in the attenuation lengths for protons 
(110 Mpc) and photons (2.6 Mpc and 2.2 Mpc [0]). At higher energy the attenuation 



length for protons dramatically decreases (13.4 Mpc aX E = 1 ■ 10^^ GeV) and the fluxes of 
protons and photons become comparable. 

A requirement for the models explaining the observed UHE events is that the distance 
between sources must be smaller than the attenuation length. Otherwise the flux at the 
corresponding energy would be exponentially suppressed. This imposes a severe constraint 
on the possible sources. For example, in the case of protons with energy E ~ (2— 3) -10^^ GeV 
the proton attenuation length is 19 Mpc . If protons propagate rectilinearly, there should 
be several sources inside this radius; otherwise all particles would arrive from the same 
direction. If particles are strongly deflected in extragalactic magnetic fields, the distance 
to the source should be even smaller. Therefore, the sources of the observed events at the 
highest energy must be at a distance i? ^ 15 Mpc in the case or protons. 

In our model the distance between sources, given by Eq.(|^), satisfies this condition for 
r > 3-10^, while Vmax ~ Vm/Vs can be many orders of magnitude larger. This is in contrast to 
other potential sources, including supeconducting cosmic strings and powerful astronomical 
sources such as AGN, for which this condition imposes severe restrictions. 

The difficulty is even more pronounced in the case of UHE photons. These particles prop- 
agate rectilinearly and their absorption length is shorter: 2 — 4 Mpc at -E ~ 3 ■ 10^^ GeV. It 
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is rather unrealistic to expect several powerful astronomical sources at such short distances. 
This condition is very restrictive for topological defects as well. The necklace model we 
introduced here is rather exceptional. 

In conclusion, we do not claim that we found a model explaining the observations, but 
almost. 
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FIG. 1. Predicted proton and gamma-ray fluxes from necklaces. The data points are fluxes from 
the compilation made in Ref.[23]. 



10 



